Coherent p° photo-production in bulk matter at high energies 
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The momentum transfer Afe required for a photon to scatter from a target and emerge as a p° 
decreases as the photon energy k rises. For k > 3 x 10 14 eV, Afc is small enough that the interaction 
cannot be localized to a single nucleus. At still higher energies, photons may coherently scatter elas- 
tically from bulk matter and emerge as a p°, in a manner akin to kaon regeneration. Constructive 
interference from the different nuclei coherently raises the cross section and the interaction prob- 
ability rises linearly with energy. At energies above 10 23 eV, coherent conversion is the dominant 
process; photons interact predominantly as p°. We compute the coherent scattering probabilities 
in slabs of lead, water and rock, and discuss the implications of the increased hadronic interaction 
probabilities for photons on ultra-high energy shower development. 

PACS numbers: PACS Numbers:13.60.-r, 25.20.-x, 95.55.Vj 
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An understanding of the cross sections for photon in- 
teractions is important in many areas of physics. For 
example, several groups have searched for radio waves 
[l[ or acoustic pulses @ produced by interactions of as- 
trophysical v e with energies up to 10 25 eV. The radio and 
acoustic frequency spectra and angular distribution de- 
pend on the distribution of moving electric charges (for 
radio) and energy deposition (for acoustic). These dis- 
tributions are controlled by the particle interactions that 
govern shower development. The experimental flux limits 
depend on the assumed character of these interactions. 

In this letter, we discuss hadronic interactions of pho- 
tons 3, and introduce a new effect: coherent photon to 
p° conversion in bulk matter. These effects become im- 
portant at energies above 10 20 eV, but are not in many 
current calculations Q. With these effects, photons 
produce hadronic showers, rather than electromagnetic 
showers. Since hadronic showers are much less affected by 
the Landau- Pomeranchuk-Migdal (LPM) effect @, they 
are much more compact; the presence of hadronic inter- 
actions alters both the frequency spectrum and angular 
distribution of the electromagnetic and acoustic radiation 
from neutrino induced showers. Published limits that do 
not consider these effects may need to be revised. 

Besides it's importance for v searches, coherent con- 
version is very interesting in it's own right, as one of a 
handful of examples where particles interact very differ- 
ently in bulk matter than with isolated atoms; interac- 
tions with individual targets are replaced by distributed 
interactions which are de-localized over multiple atoms 
Q . The other prominent examples are LPM suppression 
of bremsstrahlung and pair production 0, 0], kaon re- 
generation , and coherent neutrino forward scattering 
[111 ] . The LPM effect decreases the cross section for pair 
production, due to destructive interference. In contrast, 
in the other examples, including coherent p conversion, 



constructive interference raises the cross section. 

Photon to p conversion occurs when a photon fluctu- 
ates to a virtual qq pair which then scatters elastically 
from a target nucleus, emerging as a vector meson Q. As 
the photon energy k rises, the required momentum trans- 
fer Afc = A/ 2 7r /2fc decreases, and the coherence length 
If = h/Ak rises (we take H = c = 1 throughout). Here 
M^n is the final state mass. When k > 3 x 10 14 eV (for 
Mtttt = Mp = 778 MeV/c 2 , the p pole mass), If > 0.2 
nm, the typical internuclear spacing, and coherence over 
multiple nuclei becomes possible. 

We consider a high-energy plane-wave photon travel- 
ling in the +z direction. The photon wave function is 
mostly bare photon >), but with a qq component, 
\qq > and heavier fluctuations (e.g. qqqq, and strange 
quark pairs): 



\i/> >= e tkz [y/l-F 2 ^ > +F\qq > 



(1) 



where F = 6.0210 -2 is the amplitude for the photon to 
fluctuate to a qq pair [tJ. Heavier fluctuations are not 
relevant here, and will not be further considered. 

The \qq > scatters from nuclei at positions r*j, with 
scattering amplitude f(9), emerging with momentum k' . 
Neglecting photon or p absorption, and photon scatter- 
ing (so qq do not scatter and then fluctuate back to a 
photon), the wavefunction at a depth L in the material 
is 



^P(L) >= Vl-F^e lkL \ lb > e lkz *f(0) 
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The photon can be absorbed in the target, with absorp- 
tion length a = l/n[a ee (jA) + a ha dr('yA)]. Here a ee (jA) 
is the pair production cross section, <Thadr{lA) is the pho- 
tonuclear cross section, and n is the target density, in 
atoms /volume . 
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The p can be lost by direct hadronic interaction, or 
by decay to 7r + 7r~, followed by tt interactions. The de- 
cay itself does not affect the coherence [1^ . but two 7r 
interact differently from one p. For simplicity, we take 
(3 = l/natot{pA) where <Jtot{pA) is the p nucleus cross 
section. The factor of 2 difference between two tt and 
one p is small compared to the other uncertainties. This 
also applies for direct 7r + 7r _ production, discussed below. 

Including these absorption factors, and treating the 
target as a homogenous medium with constant density 
n, the wave function at depth L in the slab is 



< p° | ip(L) >= nF / r\dr\ I dip 
Jo Jo 
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for an incident plane wave moving along the +z axis. The 
integral over the target volume is in cylindrical coordi- 
nates where rj is the radial distance from the z axis and 
ij) is the azimuthal angle. We neglect scattering effects at 
large distances since absorption limits the effective size 
of the target. 

In an infinite medium only forward scattering, /(0), 
contributes to coherent scattering [ToL [TTL |l3j |. We as- 
sume that the real part of /(0) is small [7|, corresponding 
to nearly complete absorption. Then, from the optical 
theorem, |/(0)| = <r tot (pA)k / Air . 

After scattering, the system has the momentum k' — 
k — Ak. Since Ak cx M 2 ^, the p width affects the cal- 
culation. The photon may also fluctuate directly into 
an 7r + 7r~ pair [l4j ; these two processes interfere and the 
mass spectrum is a Breit-Wigner spectrum. Then, 
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where A and B are the energy-dependent amplitudes for 
p and direct tttt production respectively [la ], and the 
p° width T p = 150 MeV. We assume that B/A is inde- 
pendent of photon energy and target material. Integrat- 
ing J dM nw f{0,M^) returns the traditional /(0). With 
this, u{jp — * pp) = A 2 . 

We substitute r\dr\ — r'dr 1 where r' is the distance 
between the scattering point (77, z) and the observer at 
(0,-L), to evaluate the integrals (llj : 

l-M p +5r p 

< p°\iP(L) >= 2nnF / dAW(0, M^) 
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For simplicity, we give here the probability for a fixed 
(although the full calculations include the wide p): 

P p o{L,M^) = 4n 2 \f(0)\ 2 F 2 n 2 



J cos(AfcL) 
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Eq. ([6]) illustrates some of the features of the system. 
When AkL ^> 1 the cosine fluctuates rapidly, leading to 
incoherent p production. However, when the coherence 
condition AkL <C 1 is fullfilled, the scattering amplitudes 
add in phase, and production is coherent. 

A p inside a target is not directly observable. However, 
the way it interacts - electromagnetically, hadronically, or 
through p conversion - is indirectly observable, by study- 
ing the shower development. The probability that an 
incident photon interacts as a real p is the integrated 
probability for finding a p, multiplied by the probability 
for a p to interact hadronically in length dz, or dzj (3: 
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This equation is only properly normalized for P{L) <C 1. 
The loss of photon intensity due to the coherent p reac- 
tion is not included. For kaon regeneration, the analo- 
gous problem has been solved recursively [la ]. Here, we 
normalize the probabilities to sum to one in thick targets. 

The numerical values of P(L) depend on atot{pA) and 
the cross section for a photon to interact hadonically (as 
a p), (■yA) . They are related via 



(ThadrijA) = F 2 a to t(pA). 



(8) 



We determine these using two different methods: a 
Glauber calculation based on HERA data on jp — > pp, 
using a soft Pomeron model to extrapolate the cross sec- 
tions to higher energies, and a second calculation that 
includes generalized vector meson dominance (GVMD) 
plus a component for direct photon interactions [l7j . 

We consider three materials: water, standard rock Q 
and lead. For water, we add the amplitudes for the cross 
sections for hydrogen and oxygen. 

The Glauber calculation uses the optical theorem [l8j] 
and vector meson dominance model to link the total pp 
cross section to the differential jp cross section: 



vtot(pp) 



F 



Ay^TT da(jp — > pp) 



dt 



\t=o 



(9) 



where t is the squared momentum transfer. da/dt(jp — > 
pp comes from a fit to HERA data fl9j ]. 

The cross sections for nuclear targets are found with a 
quantum Glauber calculation [20| : 



a tot { P A) = 2 / ffll-e 



7totXf>P)T A (T) 



(10) 



where Ta(t) is the nuclear thickness function for a ma- 
terial with atomic number A, calculated from a Woods- 
Saxon nuclear density, with a skin thickness of 0.53 fm 
and density p = 1.16A 1 / 3 (1 - 1.16^ 2 / 3 ). 

The second approach follows Engel, Ranft and Roessler 
(ERR) [13 . It uses GVMD plus direct photon-quark 
interactions, to determine (Jhadr (jA) . ERR predict a 
steeper rise in the cross section than the Glauber cal- 
culation. We parameterize their results. For W 1P < 10 11 
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FIG. 1: The cross sections for photon interactions to e + e~ 
pairs, and for incoherent photonuclear interactions in the 
Glauber and ERR models, for water (top), standard rock 
(middle) and lead (bottom). 




10 10 2 10" 
Slab length in m 



FIG. 2: p probability as a function of depth in a slab for a 
10 23 eV photon incident on a water target. The two broad 
peaks around 20 and 40 cm correspond to AkL = 2nn at the 
p pole mass. 
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FIG. 3: The probability of a 10 24 eV photon interacting elec- 
tromagnetically (dashed blue line), in an incoherent hadronic 
interaction (dotted line), or as a coherently produced p (solid 
line), as a function of target depth in a lead target. 



eV, the cross section is constant, while at higher ener- 
gies it rises as W®p . The cross section scales as A - 887 , 
normalized so a(jPb) — 15 mb at low energies. Again, 
atot(pA) — (Jhadr{lA)/F 2 . The ERR cross sections are 
similar to a newer result that computed photon cross sec- 
tions using a dipole model [2lJ . 

Figure [T] compares the cross sections for e + e _ produc- 
tion and the two photonuclear models. a ee {^A) is con- 
stant at low energies, but falls at higher energies when the 
LPM effect becomes important [!, H[ . The hadronic cross 
sections rise slowly with energy. The hadronic curves 
agree at low energies, but diverge as the energy rises. In 
the ERR (Glauber) model, above 5 x 10 19 eV (5 x 10 20 
eV), photonuclear interactions predominate, and photons 
produce hadronic showers. The cross-over energy is al- 
most material-independent for solids. 

Figure [2 shows the probability of finding a p°, Eq. (J6]), 
as a function of depth in a thick target. As with K s 
from kaon regeneration [Io| , the amplitudes add and the 
probability initially rises as the square of the depth. After 
reaching a maximum, the probability decreases slowly as 
p absorption competes with p production, and then drops 
rapidly as the bare photons are absorbed. 

Figure [31 shows the probability for a 10 24 eV photon to 



interact as a function of target thickness L. For L < 10 
cm, the probability of a photon interacting as a coher- 
ently produced p is quite small. For L > 10 cm, the 
probability is higher than for photonuclear interactions 
and e + e~; it is the dominant interaction! 

Figure 0] shows the probabilities for incident photons 
to interact electromagnetically (7 — * e + e~ pair), via in- 
coherent hadronic interactions, and as coherently pro- 
duced p, as a function of energy, in a 100 m thick slab; in 
the displayed energy range, 100 m is effectively infinite 
thickness. At energies above 10 23 eV (for ERR) to 10 24 
eV (for the Glauber calculation), coherently produced p 
interactions dominate. This dominance will reduce the 
interaction length of photons. The reduction in shower 
length will have consequences for searches for ultra-high 
energy astrophysical u e |l|, |4[ . 

In conclusion, photons with energies above 10 20 eV in 
solids or liquids are more likely to interact hadronically 
rather than electromagnetically. These hadronic showers 
are not subject to the LPM effect, and so are considerably 
more compact than purely electromagnetic calculations 
would predict. At energies above 10 23 eV, photons are 
most likely to interact by coherently converting into a p, 
and then interacting hadronically in the target. This is a 
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FIG. 4: The normalized probabilities for a photon to interact as an e + e~ pair, via incoherent hadronic interaction ('hadron'), 
or as an elastically produced p ('rho') in a 100 m thick water or standard rock target. This is thick enough for almost total 
absorption, so these results apply for an infinitely thick target. 



new example of a coherent process in bulk matter, simi- 
lar to kaon regeneration or coherent neutrino scattering. 
This coherent interaction further shortens the shower de- 
velopment. The reduction in photon interaction length 
will shorten ultra-high energy electromagnetic showers, 
altering the radio and acoustic emission. It should be 



included in the models used in searches for ultra-high 
energy astrophysical u e . 
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